One-dimensional ballistic transport with FLAPW Wannier functions 
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We present an implementation of the ballistic Landauer-Biittiker transport scheme in one- 
dimensional systems based on density functional theory (DFT) calculations within the full-potential 
linearized augmented plane-wave (FLAPW) method. In order to calculate the conductance within 
the Green's function method we map the electronic structure from the extended states of the FLAPW 
calculation to Wannier functions which constitute a minimal localized basis set. Our approach ben- 
efits from the high accuracy of the underlying FLAPW calculations allowing us to address the 
complex interplay of structure, magnetism, and spin-orbit coupling and is ideally suited to study 
spin-dependent electronic transport in one-dimensional magnetic nanostructures. To illustrate our 
approach we study ballistic electron transport in non-magnetic Pt monowires with a single stretched 
bond including spin-orbit coupling, and in ferromagnetic Co monowires with different collinear mag- 
netic alignment of the electrodes with the purpose of analysing the magnetoresistance when going 
from tunneling to the contact regime. We further investigate spin-orbit scattering due to an impu- 
rity atom. We consider two configurations: a Co atom in a Pt monowire and vice versa. In both 
cases, the spin-orbit induced band mixing leads to a change of the conductance upon switching the 
magnetization direction from along the chain axis to perpendicular to it. The main contribution 
stems from ballistic spin-scattering for the magnetic Co impurity in the non-magnetic Pt monowire 
and for the Pt scatterer in the magnetic Co monowire from the band formed from states with d xy 
and d x 2_ y 2 orbital symmetry. We quantify this effect by calculating the ballistic anisotropic mag- 
netoresistance which displays values up to as much as 7% for ballistic spin-scattering and gigantic 
values of around 100% for the Pt impurity in the Co wire. In addition we show that the presence 
of a scatterer can reduce as well as increase the ballistic anisotropic magnetoresistance. 

PACS numbers: 



I. INTRODUCTION 



With the possibility to perform transport measure- 
ments on nano- down to atomic-scale junctions using 
mechanically-controllable break-junctions^ or scanning 
tunneling microscopjJ^El various fundamental questions 
on electron transport as well as practical problems con- 
cerning device functionality have arisen. With shrink- 
ing system size the junctions have become considerably 
smaller than the mean free path of a transmitted electron, 
reaching the ballistic transport regime. In this regime 
various effects such as the geometric arrangement of the 
atoms, the chemical composition, the magnetic order, vi- 
brations, correlation effects, or the magnetic anisotropy 
can play an important role due to the reduced coordina- 
tion number of the participating atoms. In the context of 
spin-dependent transport, for example, there is a strong 
interest in understanding how the spin- valv e effect scales 
to systems of atomic- or molecular-scale^'. In nanoscale 
junctions, new transport effects can also arise suc h as 
the ballistic anisotropic magneto-resistance (BAMRp^. 
In order to successfully address such issues a theoretical 
description needs to properly take into account the elec- 
tronic structure of the system which is typically obtained 
by first-principles methods based on density functional 
theory. The central experimental quantity is the mea- 




FIG. 1: (color online) Typical geometry of a ballistic trans- 
port calculation, consisting of three different regions (left 
lead/scattering region/right lead). Charge carriers with en- 
ergy E are transmitted through a scattering region with a 
transmission probability T(E) from the left lead to the right 
lead. The blue planes separate the leads from the scattering 
region. While the semi-infinite leads resemble the electronic 
structure of a periodic system, the scattering region includes 
the scatterer as well as the lead-scatterer contact region. 

sured current versus bias voltage (I-V-curve) or at small 
bias voltages the conductance. 

The theoretical method most frequently applied to 
describe quantum transport in such systems is the 
Landauer-Biittiker approach in which the junction is di- 
vided into a central scattering region and two leads in 
thermal equilibrium with contact reservoirs, Fig. [TJ re- 
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suiting in the famous Landauer conductance formulc^. 
Basically two different groups of techniques have been 
developed to solve the transport problem: wave-function 
and Green's function (GF) based methods which are 
equivalent in case of non-interacting charge carriers^. 
Among the wave function based methods the trans- 
mission through such a quasi one-dimensional (ID) 
system can be calculated by means of the transfer 
matrix method^^, solving the Lippmann-Schwinger 
equatior l 18 * 19 !, or by wave function matching^. The GF 
methods are usually based on Keldysh, Kadanof f and 
Baym's non-equilibrium Green's functions (NEGFjPMH. 
Beyond the standard non-interacting electron approach, 
there has bee n wor k incorporating e.g. inelastic scattering 
on vibration d 23 ! 24 ! or treating correlation effects through 
self-energies 25 . An alternative way of calculating quan- 
tum transport is by using the Kubo approach as formu- 
lated by Baranger and St one^ , relating the current to 
the dynamical polarizatiorPZGS. 

Based on these three general approaches all codes differ 
in the way the electronic structure is described. In the 
first implementations based on density functional theory 
(DFT), the electrodes were tre ated as j ellium which were 
coupled to the scattering re gi or ]1511Z118| Large systems up 
to devices can be described using semi-em pirical tight- 
binding methods for the electronic structur d 13 * 14 * 19 * 2 ^^ 
while approaches using DFT for both the description of 
the electrodes via self-energi es and t he scattering region 
promise the highest accuracjP^SEEffll 

Among these implementations various DFT meth- 
ods have been applied. Transport codes based on 
Green's functions rely on a localized basis set limit- 
ing this approac h to basis sets of nu merical orbit als 
such as Gaussians^Hl, localized orbit alfpEMHIIII] or 
wavelets^. The application of flexible and accurate 
plane-wave DFT methods for transport calculations is 
usually realized in connect ion w ith the scattering ap- 
proach for the conductanc d 46 * 47 !. Alternatively, the ef- 
ficient GF method for the transport calculation can 
be used if the extended states in the plane-wave ex- 
pansion are mapped onto maximally localized Wannier 
functions^ (MLWFs). This approach combines plane- 
wave calculations with the use of a mini mal basis set 
suitable for quantum transport calculation a 36 ! 37 " 3 " 9 ! 50 ! . 
When one is dealing with low-dimensional systems and 
subtle band structure effects such as spin-orbit coupling 
the accuracy of electronic structure description becomes 
crucial. Therefore, the application of a highly precise all- 
electron full-potential linearized augmented plane-wave 
code is desirable. To our knowledge no such DFT trans- 
port scheme has been reported an d only a few codes allow 
to incorporate spin-orbit coupling} 47 * 51 * ^. 

In this paper we present a method to calculate 
transport through ID nanoscale structures following 
the Landauer-Buttiker approach. The underlying elec- 
tronic structure of the studied system is obtained from 
DFT within the ID version^! of the full-potential lin- 
earized augmented plane- wave method (FLAPW), as im- 



plemented within the FLEUR codeP The 1D-FLAPW 
method is specifically tailored to treat ID structures 
avoiding supercell calculations: the periodicity is explic- 
itly taken into account only along the nanostructure's 
axis (z-axis in the following), while the wavefunctions in 
the vacuum surrounding the system are forced to obey an 
exponential decay in-planePS Since the FLAPW wave- 
functions are intrinsically delocalized in real space we 
perform a mapping of the electronic structure of the sys- 
tem onto a set of localized Wannier functions (WFs), 
which allows to solve the transport problem in real space 
efficiently. The WFs obtain ed from the FLAPW calcu- 
lation (FLAPW WFsp 515 ^ provide a minimal localized 
basis set which describes the ab initio electronic bands 
within a certain energy window with high accuracy and 
allows to efficiently compute the non-equilibrium Green's 
function (NEGF) of the system needed to determine its 
transmission function T(E). We use and compare two 
different sets of WFs, namely, the maximally-localized 
Wannier functions (MLWFs), 48 which are uniquely de- 
fined by fulfilling the condition of maximal localization in 
real space, and the so-called first-shot Wannier functions 
(FSWFs)P^ being much easier to obtain computation- 
ally and although non-unique, still capable of describing 
the transport properties of a system correctly in many 
cases. A special approximation we include in our trans- 
port scheme is the so-called "locking technique", which 
allows to use separately calculated leads and scattering 
regions and to combine those into one quantum transport 
calculation, achieving an accurate treatment of leads and 
scattering region at reduced computational cost. 

As a first application we have calculated the electronic 
structure and the ballistic transport properties of a non- 
magnetic Pt monowire with a single stretched bond in 
the middle of the chain which acts as a source of scatter- 
ing. For this rather simple system we demonstrate the 
quality of our MLWFs and FSWFs, the locking technique 
to obtain the Hamiltonian of the open system, and show 
the possibility of decomposing the transmission function 
in terms of orbital symmetry. We further investigate the 
influence of spin-orbit-coupling (SOC) on the transmis- 
sion of the Pt wire. We find a substantial change of the 
conductance of one quantum of conductance at the Fermi 
level for a perfect wire due to the strength of SOC in 5d- 
transition metals such as Pt. 

In order to include the effect of large spin-polarization 
we have chosen a ferromagnetic Co monowire with a 
single stretched bond, a prototypical magnetic system, 
and calculate the magnetoresistance from the conduc- 
tance in a parallel and antiparallel alignment of the Co 
electrodes. We obtain a rapid decrease of the mag- 
netoresistance with the separation between the two Co 
monowires which is due to the fast decay of transmission 
from the highly spin-polarized localized states of d xz ^ yz — 
and d X y tX 2_ y 2— symmetry. 

Finally, we have studied scattering from a single im- 
purity atom in a monowire due to SOC. We have chosen 
two configurations: (i) a non-magnetic Pt atom in a fer- 
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romagnetic Co wire and (ii) a magnetic Co atom in a 
non-magnetic Pt wire. In both cases, we have compared 
the conductance obtained in the scalar-relativistic ap- 
proximation and upon including SOC. We find a strong 
influence of SOC on the transmission due to the induced 
splitting of bands. In addition, the conductance depends 
sensitively on the magnetization direction in the system 
being either along the wire axis or perpendicular to it. 
While in case (ii) the resulting ballistic anisotropic mag- 
netoresistance displays values of 7 % due to spin-orbit in- 
teraction mediated scattering into both spin-channels for 
the symmetry breaking out-of-chain quantization axis, 
in case (i) the values of BAMR reach as much as 100%, 
reflecting the giant value of the ballistic anisotropic mag- 
netoresistance of the pure Co chain. 

The paper is organized as follows. In Sec. [TT] we de- 
scribe the theoretical basis of our approach to calculate 
the conductance and introduce the key quantities. In 
particular, the Green's function method is applied to ob- 
tain the transmission function and the conductance. The 
mapping of the electronic structure from the FLAPW 
method to a localized basis set is accomplished via Wan- 
nier functions. The construction of the Hamiltonian for 
the open quantum system is described. In Sec. 



where H 



L/R 



is the semi-infinite Hamiltonian of the 
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present the first applications of our new transport code 
to several typical systems of interest. We begin with the 
conductance for a non-magnetic Pt wire with a single 
broken bond and study the transmission as a function 
of bond length and upon including spin-orbit coupling. 
Then the magnetoresistance of Co monowires with a sin- 
gle elongated bond is discussed. Finally, the effect of 
spin-orbit scattering is illustrated by two examples: a 
Pt monowire with a single magnetic Co atom and a Co 
monowire with a single Pt atom. A summary and con- 



clusions are given in Sec. VI 



II. METHOD 
A. General transport problem 

We describe the transport properties of the system 
within the Landauer-Buttiker approach, dividing it into 
three different regions: two semi-infinite leads (left, L, 
and right, R) and the scattering region (S), which in- 
cludes the actual scatterer as well as the lead-scatterer 
contact region, in which the effect of the scatterer on 
the properties of the leads ideally decays such that their 
electronic structure can be considered perfect and unper- 
turbed inside the L and R regions. Assuming that the 
interaction between the left and right leads can be ne- 
glected the tight-binding Hamiltonian of our system cor- 
responding to such a structural division has the following 
form: 
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left/right lead, while Hls/SR describes the coupling of 
the scattering region to the leads and H5 is the Hamilto- 
nian of the scattering region. Due to the semi-infiniteness 
of the leads the dimension of the Hamiltonian, Eq. ([I]) , is 
infinite, which presents a conceptual computational prob- 
lem. An efficient method to deal with that, applicable to 
any system of the type depicted in Fig. 1, which can be 
described with a real-space tight-binding Hamiltonian of 
the type of Eq. (1), has been developed.^ This method is 
based on the non-equilibrium Green's function (NEGF) 
formalism, which treats the scattering region and the 
semi-infinite leads on equal footing, describes extractions, 
re-injections and excitations of electrons in the system 
and solves the problem of the semi-infinite leads by in- 
troducing finite-dimensional self-energies H L / R which in- 
clude the true lead's effect on the scattering process. 
Within the NEGF formalism the system is described by 
means of the retarded Green's function: 



G(E) = [{E- 



)I H1 



(2) 



where I denotes the unity matrix of the dimension of H. 
By neglecting at first the coupling of the leads to the scat- 
tering region and regarding just the first few layers of the 
leads which are actually interacting with the scattering 
region, it is possible to replace the leads Green's function 
by their surface Green's functions §,l/r{E)W^ This can 
be derived by rewriting the lead's Hamiltonian in a block 
diagonal form using square matrices h. L / R and h. LL / RR 
of the same dimension as the surface Green's function: 
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( ■. \ 

h LL h L h[ L 
V h LL h L ) 



(3) 



Based on this description of the leads, the surface 
Green's function Sl/r(E) can be determined iteratively 
starting from 



gf /R (E) = [(E + ie)I L/R 



h L/R ]-\ 



(4) 



with g^ R (E) being a square matrix with the dimen- 
sion of interacting orbitals at the leads' surfaces. The 
expression Q can be converged to the surface Green's 
function by recursively incorporating the inter-layer in- 
teraction sub-matrices h LL i RR with an efficient recursive 
scheme P^l 

By reintroducing the coupling of the scattering region 
to the leads as a perturbation to the system, the Green's 
function Gs(E) of the scattering region can be obtained 
from the unperturbed Green's function of the scattering 
region by the Dyson equation 



G S (E) = [Els - Us - H T LS g L H LS - H^gflHsij]" 1 . 

(5) 

The whole effect of the semi-infinite leads on the con- 
ductor can be then expressed by the leads' self-energies 
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Hl/r(E), which incorporate the surface Green's func- 
tion gi/ij(-B) and the now finite-sized coupling matrices 
H-LS/SR) adapted to the size of the surface Green's func- 
tions: 

^l/r{E) = H^ LS/SR g L/R (E)B. LS/SR . (6) 

The self-energies are obviously finite-sized matrices of the 
dimension of H5. The self-energies are related to the 
broadening matrices T: 

T L/R (E) = i{V L/R {E) - Vl /R (E)}, (7) 

which describe the effect of broadening of the states in 
the scattering region caused by the presence of the leads 
as well as the transfer rates of charge carriers from the 
leads into the scattering region. The incorporation of the 
non-Hcrmitian self-energies changes the nature of the de- 
scription from the static steady-state picture of the open 
system to a dynamic transport scheme, responding to an 
incoming charge carrier with the energy E. Based on 
these quantities the transmission function T(E) , describ- 
ing the probability of charge carriers originating from one 
lead to be transmitted to the other lead, can be expressed 
in the following way: 

T(E) = TiiG s (E)T L (E)Gl(E)T R (E)}. (8) 

The current, being a natural observable in a quantum 
transport measurement, can then be calculated from the 
Landauer formula: 

I= P -J dET(E)[f L (E) - f R (E)j, (9) 

where Jl/r are the occupation functions of the leads. 
The expression for the conductance then reads: 

G(E) = ^T(E). (10) 

In case of perfect transmission, T(E) — 1, this results in 
the well-known conductance quantum 

2e 2 

Ob- x . (ID 
for a single, spin-degenerate band. 

B. Prom FLAPW states to localized Wannier 
functions 

The aim of the approach introduced here is to combine 
the accuracy and speed of state-of-the-art DFT electronic 
structure calculations based on the one-dimensional ver- 
sion of the FLAPW method as implemented in the FLEUR 
codej^ and the capability of the NEGF-formalism de- 
scribed above to treat the whole transport problem in an 
efficient way. Especially for transport phenomena driven 



by magnetism or spin-orbit coupling (SOC) a precise de- 
scription of the electronic structure is necessary. Typical 
systems currently under scrutiny in experiment include 
geometries with a low coordination number which favors 
magnetism and gives rise to strong SOC due to unquench- 
ing of the orbital momentP 

The major problem in combining an LAPW or a plane- 
wave based electronic structure method with the real- 
space transport schemes lies in the fact that normally 
several hundreds of delocalized basis functions per atom 
are used in such codes in order to achieve the required 
accuracy. In our implementation we use the machinery 
of Wannier functions (WFs) , constructed out of FLAPW 
wavefunctions J^l which proved to be an efficient connec- 
tion between the two, conceptually independent, compu- 
tational methods. The main advantage in such a "link" 
can be attributed to two factors: (i) using the gauge 
freedom of Wannier functions they can be enforced to be 
rather localized in real space, and (ii) an "exact" map- 
ping of the ab initio Hamiltonian onto a tight-binding 
representation with WFs as a localized orthonormal ba- 
sis set can be achieved.^ 

Having at hand the converged Bloch wavefunctions 
■0mk for a set of bands m < M calculated on a uniform 
mesh of N fc-points, the orthonormal set of Wannier func- 
tions can be obtained via the following transformation^ 

M 

\w Rn ) = jfY, e ~ ik R E u ™ ivw) , (12) 

k m— 1 

where the number of WFs N should be smaller than or 
equal to M. The gauge freedom of WFs manifests itself 
in that the matrices U^ ln (in the following, U-matrices) 
can in principle be arbitrary. In the case when N = M 
and the group of bands from which we are extracting the 
WFs from is isolated from other bands, the U-matrices 
are unitary at each fc-point. Imposing the constraint of 
maximal localization of WFs in real space determines 
the set of U-matrices up to a common global phase, and 
the corresponding set of WFs is called the maximally- 
localized Wannier functions (MLWFs) For the whole 
procedure of maximal localization we use the Wannier90 
codeP 

The criterion for the localization of WFs is the small- 
ness of their spread!^ The process of the spread min- 
imization constitutes an iterative process at the end of 
which the U-matrices corresponding to the MLWFs are 
obtained. This minimization procedure requires as a 
starting point a certain initial guess for the set of the 
MLWFs. In order to construct this set, one chooses cer- 
tain localized orbitals \g n ), which are projected onto the 
subspace of wavefunctions \ifrk m ): 

m 

and then orthonormalized: 

^k„)=£(( S(k) )pl^ m >, (14) 

m 
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with the overlap matrix Sm„ = (</>kn 
the starting WFs can be generated: 



k 



V'kn 



after which 



(15) 



This orthonormal set of Wannier orbitals we will call in 
the following the first-shot WFs (FSWFs). 

The FSWFs are not unique in the sense that they 
strongly depend on the choice of the localized orbitals 
g n . In many cases, however, especially when MLWFs 
are well-localized around atoms as in the case of certain 
rf-orbitals in most of transition-metals and transition- 
metal oxides,^ the difference between the FSWFs, orig- 
inated from the localized d-orbitals, and the correspond- 
ing MLWFs is rather small. This allows to spare the 
computational time needed for the minimization of the 
spread, and immediately construct, e.g., the needed effec- 
tive Hamiltonians in terms of FSWFs. Examples, when 
there is a substantial difference between the FSWFs and 
MLWFs, include orbitals for which the centers of the 
WFs do not coincide with the centers of atoms. In the fol- 
lowing we will analyze in detail the difference in transport 
properties calculated with MLWFs and FSWFs, both in 
the case when there is little difference between the two 
sets of WFs and when the difference between them is 
significant. 



C. Construction of the Hamiltonian in real space 

In terms of the FLAPW basis functions the Hamilto- 
nian can be written as 



HflAPW = X^m(k) \lpmk) (V>rok| 



(16) 



while in terms of WFs the equivalent expression is 

H WFs = J2J2 #™,«'( R 1 - R 2) IWnlO (Wn'Rj , 



where 

fln,n'(Rl ~ R 2 ) = (W nRl | H FL APW I WVr 2 ) (18) 



are the hopping integrals between the n and n' Wannier 
orbitals at sites R4 and R2. By substituting Eq. (16) 
into Eq. ( 18 1 we find: 



(R4-R2) = ~ ^£ m (k) (W„ Rl I -0mk) (V>mk I W n > R2 

m,k 

(19) 

Thus the real-space representation of the Hamiltonian 
in terms of localized Wannier functions can be derived 
from the knowledge of the eigenvalues and wavefunctions 
of the system. In respect to WFs, for efficient evalua- 
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FIG. 2: (color online) Schematic sketch of a ballistic trans- 
port calculation based on a WFs tight-binding Hamiltonian. 
The leads are described by perfect wires to exclude spuri- 
ous deviations from their exact electronic structure. Their 
semi-infinite structure is constructed from the Hamiltonians 
of principal layers h.L/R and the interaction matrices h. LL / RR 
between two principal layers. The scattering region is de- 
scribed by the Hamiltonian Hj and coupled to the leads by 
the interaction matrices ~H LR /sr, extracted from a supercell 
calculation. The supercell has to be large enough to repro- 
duce the lead-scatterer contact with desired accuracy, usually 
larger than sketched here. 

required!^ The correspondence in the eigenspectrum be- 
tween the constructed Hamiltonian in terms of WFs and 
the Hamiltonian in terms of eigenfunctions is exact on 
the grid of fc-points used for the WFs construction, and 
for this reason the set of WFs is sometimes referred to as 
exact basis set, or, the tight-binding basis set of ab ini- 
tio accuracy. (This is only valid within the frozen inner 
window of disentangled systems^). 

According to Eqns. Q to (8J), in order to calculate the 



(17) transmission function only t 



re h 



L/R- 



ilLL/RR, H-LS/SR 



and H5 parts of the Hamiltonian are needed. Given a 
FLEUR electronic structure calculation it is necessary to 
construct these parts of the Hamiltonian from the result- 
ing WFs hopping elements (Eqn. 
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We focus on the 

correct treatment of the scattering region (see Fig. ([2])). 

After determination of the atoms belonging to the scat- 
tering region, it is possible to write down the preliminary 
result for H5, based on Eqn. (19 1, 



(20) 



tion of Eqn. (191 only knowledge of the U-matrices is 



where i and j determine the atom and n and m the in- 
herent WFs. 

Due to the real-space decay of the WFs the correspond- 
ing hopping matrix elements i7„.„/(Ri — Rj) also decay 



(> 



as the distance in real space between the Wannier func- 
tions |Rj — Rj | is increasing. For an efficient use of the 
real-space WFs Hamiltonian within the transport scheme 
described above it is necessary to keep its matrix ele- 
ments only up to a certain number of nearest neighbors 
(n.n.), setting the rest of the elements to zero. As a re- 
sult of this procedure the Hamiltonian matrix becomes 
sparse, which allows for a computationally inexpensive 
treatment. For a given number of n.n., the quality of 
the sparse Hamiltonian depends on the degree of local- 
ization of the WFs. Here, by quality of the Hamiltonian 
we mean the correspondence between its eigenvalue spec- 
trum to that obtained from ab initio, or, in the sense of 
transport, how well-converged the transmission function 
T(E) is with respect to the number of n.n. In this respect, 
in the following we compare and analyze the results ob- 
tained with MLWFs and FSWFs, which display different 
localization properties. 

One way to deal with the exponential decay in 



Eqn. (20 1 would be to manually eliminate all matrix ele- 



ments beyond a certain n.n.. We propose here a flexible 
scheme, minimizing this effort by dividing the scattering 
region into principal layers h/, I = 1, . . . , s and interac- 
tion matrices h^+i between neighboring layers: 



He 



hi2 
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(21) 



The sub-matrices are set up as Eqn. (20). For the on-site 
matrices, h;, the indices i and j are restricted to atoms 
from the given layer I. For the interaction matrices hu+i 
the index i is restricted to atoms from layer I and the 
index j to atoms of the neigboring layer I + 1.. 

While still capable of describing the system in terms 
of Eqn. (20) (with s = 1), the principal layers can opti- 
mally contain the number of atoms effectively interact- 
ing, reducing the number of neglected hoppings. Typi- 
cally these principal layers are chosen to contain the same 
number of atoms as the principal layers of the leads re- 
sulting in the same aproximation in terms of n.n. for both 
regions and thereby avoids inconsistencies in the transort 
calculations. Furthermore this scheme allows possible fu- 
ture extensions, such as e.g. a combination of separately 
calculated scatterers into one scattering region. 

Knowing the Hamiltonian H5 of the scattering region, 
it is necessary to determine the coupling of the scattering 
region to the leads. Since the Hamiltonians of both leads 
and the scattering region are partitioned into principal 
layers (see Eqn. (|3| and Eqn. (21)), we only need to find 
the hopping elements between the adjacent layers. Inter- 
actions between non-neighboring layers are neglected by 
construction. The non-zero elements of H^s can now be 
extracted from the supercell calculation (see Fig. [2]) as 
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(22) 



where the index i runs over the atoms of the principal 
layer 1 of H5 and the index j runs over the principal layer 
of the left lead. Hsr can be constructed analogously. To 
prevent a significant systematical error, it is necessary 
to make the original supercell large enough to screen an 
unphysical inter-unit cell interaction. 

Finally, only the Hamiltonians for the leads are miss- 
ing. Ideally, the calculated unit cell should be large 
enough in order to reproduce the properties of the bulk 
material far away from the scatterer and thus, the lead 
Hamiltonian can be extracted directly from the supercell 
calculations in a straight-forward manner. Owing to the 
significant computational burden, it is, however, hardly 
feasible to apply this approach to large and complex sys- 
tems while keeping at the same time the accuracy neces- 
sary to capture the main energy scales of the phenomena 
studied. The technique we use to overcome this problem, 
particularly prominent for the FLAPW method with its 
complicated basis set, is discussed in the following sec- 
tion. 

Up to now no comments have been made concerning 
the way magnetic systems and the effect of SOC are 
treated. For magnetic systems the majority and minority 
spin channels can be regarded separately, resulting in two 
independent calculations of the transmission function for 
spin-up and spin-down channels. In the presence of SOC 
the whole methodology holds considering that both spin 
channels have to be treated together, thus resulting in 
twice the number of WFs used simultaneously to solve 
the transport problem. 



D. Locking technique 

The accurate treatment of the leads within the ap- 
proach described above constitutes a considerable chal- 
lenge. Taken from a self-consistent supercell electronic 
structure calculation as they are, the sub-matrices h. L / R 
and h^/^/j will contain deviations from "ideal" -lead 
matrix elements in a large vicinity of the scattering re- 
gion. While some of these deviations are definitely physi- 
cal in their origin due to a large decay length of ID charge 
perturbations caused by the scatterer, the rest of them 
will be a spurious artefact of the supercell approach owing 
to the fact that the leads as calculated are not intrinsi- 
cally semi-infinite. This presents a considerable problem 
in particular when the leads have to be described with 
Hamiltonians beyond the 1 st n.n. In this case to describe 
the semi-infinite leads precisely one would have to go to 
huge supercells so that the A atoms in the supercell de- 
scribing the lead would be exactly identical, with A being 
the number of atoms in one principal layer (see Eqn. ([3|). 
We found that condition impossible to achieve for non- 
trivial systems. Another approach of constructing a lead 
beyond 1 st n.n. artificially from the outmost atoms of the 
scattering region by periodically expanding it is flawed, 
too, due to the unknown unperturbed hopping matrix el- 
ements beyond 1 st n.n.. This is a serious problem, since 
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the lead has to be described as precisely as possible to 
prevent a huge systematic error. 

The basic idea to work around this problem is as simple 
as effective, namely matching the supercell hopping ma- 
trix elements to those of the true leads. Within this so- 
called "locking" technique, the leads are replaced by the 
perfect wires, providing correct self-energies and Fermi 
levels of the true infinite periodic system, while the su- 
percell size is chosen large enough to describe the lead- 
scatterer interface region sufficiently well, see Fig. [2] In 
our transport approach this means that different parts of 
the Hamiltonian (Eqs. p0| ) and fl22j) are extracted from 
two different DFT calculationsPW the H s and H L s/SR 
coupling matrices are taken from the supercell calcula- 
tion describing the scattering region, while the / R and 
h-LL/RR sub-matrices (needed in Eqn. (|3|)) are taken from 
the calculation for the perfect leads, iil/r an d ^ll/rr. 
can be determined similar to the principal layers h; and 
h/,; + i of Hs (Eqns. ( 20 ) and ( 21 )) with the principal layer 
I and the neighboring identical layer I + I. To achieve 
matching Fermi levels for lead and supercell calculations, 
it is additionally necessary to align the diagonal elements 
of the matrices h L / R (Eqn. [3|) and Hs (Eqn. 21). 



MLWFs 



FSWFs 



III. Pt MONOWIRES 

In the following sections we present a few instructive 
applications which illustrate the quality and possibilities 
of our FLAPW-WF based approach to obtain the con- 
ductance in one-dimensional magnetic systems within the 
Landauer coherent transport method. In this section, we 
focus on Pt monowires which possess a single stretched 
bond that acts as a scattering potential for electrons. 
Starting from the construction of the WFs and the tight- 
binding like Hamiltonian, we discuss the transmission 
function and its decomposition in eigenchannels. Our re- 
sults further demonstrate the applicability of the locking 
technique described above. Finally, we include spin-orbit 
coupling in our calculations and show that the obtained 
transmission compares well with that calculated based 
on the scattering approach in combination with a pseu- 
dopotential method for the electronic structure^. 

In order to calculate the conductance within the ap- 
proach described in the previous sections we need to 
perform two separate DFT calculations and subsequent 
Wannierizations for every system: (i) a calculation for 
the semi-infinite electrode and (ii) a supercell calculation 
which includes the scattering center. From the latter, we 
determine the hopping matrix elements for the coupling 
to the leads. For the monowires considered in the follow- 
ing the Hamiltonian of the semi-infinite electrode can be 
obtained from a calculation with one atom in the unit 
cell. For the scattering region, we have used supercells 
of different size as described in the computational details 
section in the appendix. 
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FIG. 3: (color online) Comparison of the Pt monowire band- 
structure in SR approximation calculated within DFT (big 
black dots) and obtained from the WF Hamiltonian based on 
(a) MLWFs and (b) FSWFs considering a limited number of 
nearest neighbors. In (a) the orbital character of the states 
is given and in (b) the bands are denoted according to their 
symmetry with respect to the chain geometry. 



A. Bandstructure and hoppings 

Before proceeding into the discussion of the transmis- 
sion, it is insightful to examine the localization proper- 
ties of typical MLWFs and FSWFs which we use for our 
transport calculations. While the unique MLWFs are 
rather well localized in real space, this is not necessarily 
the case for the FSWFs, which strongly depend on the 
choice of the initial orbitals. If the trial orbitals do not 
differ very much from the final result of the localization 
procedure the difference in spread between the MLWFs 
and the FSWFs can be small. 

For transition-metal monowires this is the case for the 
localized d-orbitals of A3 symmetry (d xz and d yz ) and of 
A4 symmetry {d xy and d x 2_ y 2). Taking an infinite pe- 
riodic Pt monoatomic chain with an interatomic spacing 
of 4.48 bohr as an example we calculate the spreads of 
the A 3 and A 4 MLWFs to be 3.70 bohr 2 and 2.22 bohr 2 , 
respectively. The calculated spreads of the FSWFs, con- 
structed with solutions of the radial equation for the ac- 
tual potential obtained from the first-principles calcula- 
tion^, are indistinguishable from the former. 

The situation is completely different, however, for the 
FSWFs constructed from the s- and e^-like trial or- 
bitals. In this case the difference in spread between the 
resulting FSWFs and the Ai-like MLWFs is remarkable. 
While values of 2.89 and 6.20 bohr 2 are obtained for the 
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spread of d z 2-like and s-like MLWFs, respectively, the 
corresponding values constitute 55.78 and 319.44 bohr 2 
for FSWFs. This indicates that the MLWFs differ signif- 
icantly from the trial functions. 

The reason for the rather large spreads of the FSWFs 
can be found by comparing the FSWF centers to the ML- 
WFs centers. In case of MLWFs the centers of the s-like 
WFs are located between the atoms, forming covalent 
bridge-like Wannier functions. Such Wannier functions 
are hard to construct directly from the atom-centered 
trial orbitals. The FSWFs constructed from the s- and 
d z 2-like trial orbitals are, in contrast, located on the 
atoms, which causes a significantly larger spread^. 

In principle, all Hamiltonians obtained by mapping to 
Wannier functions which include the hopping matrix ele- 
ments between all WFs are equivalent. This equivalency 
is lifted, however, if we consider only a limited number 
of neighbors to set up our tight-binding like Hamilto- 
nian. In Fig. [3j the Pt monowire bandstructure based 
on the FLAPW-DFT calculation and Slater-Koster in- 
terpolations of the bandstructure based on MLWFs and 
FSWFs are compared. The trial orbitals for the FSWFs 
are in this case chosen to be s- and d-like orbitals and 
centered on each atom. While in first n.n. approximation 
the interpolated bandstructures differ between the ML- 
WFs and FSWFs approach, especially in the bandwidth 
of the more delocalized s and d z 2 orbitals, already in sec- 
ond n.n. approximation both WFs basis sets describe the 
FLAPW-DFT bandstructure equally well. By further in- 
creasing the considered number of neighbors to third n.n. 
approximation, the accuracy of the description increases 
with respect to the s-bandwidth. However, the most im- 
portant part with respect to transport properties is the 
bandstructure in the vicinity of the Fermi level, which 
does not improve significantly. For the more localized 
d xy , d x 2_ y 2, d xz and d yz orbitals even the first n.n. de- 
scription is sufficient as seen in the bandstructure and 
also from the hopping matrix elements as seen in Fig. [4] 

At least for a rather simple system such as a perfect Pt 
monowire, the localization procedure used to obtain ML- 
WFs obviously does not influence the localized <i-orbitals 
mentioned above. Only the s and d z i states are affected, 
but the decay of the hopping integrals is exponential ir- 
respective of the description (FSWFs or MLWFs). For 
systems more complicated than a Pt monowire, the ini- 
tial choice of trial orbitals may not be straightforward. 
In such a case the localization procedure to obtain ML- 
WFs significantly improves the accuracy of the calcu- 
lation, while for simpler systems where more intuitive 
choices of orbitals can be made FSWFs may be suffi- 
cient. An example that both descriptions indeed lead to 
very similar results with respect to transport calculations 
is shown below for a Pt monowire with one elongated 
bond. Note that the FSWFs make the construction of 
the transport Hamiltonian, as discussed in section [II C[ 
much more simple, expecially for systems with a more 
complex electronic structure. 
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FIG. 4: (color online) Real space hopping integrals between 
orbitals of the same type \H(Ri — Ro)| as a function of the 
n.n. for a Pt monowire on a logarithmic scale. The hoppings 
were calculated both with MLWFs (open red symbols) and 
FSWFs (closed black symbols) for (a) s and d z 2 orbitals, and 
for (b) d xy , d x 2_ y 2, d x 
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FIG. 5: (color online) Conductance for a Pt monowire with a 
single bond stretched by A = 0.72 bohr using MLWFs within 
the nearest-neighbor approximation for the transport Hamil- 
tonian and (a) a 6 atom supercell and (b) a 12 atom super- 
cell for the FLAPW calculation of the scattering region. The 
semi-infinite leads have been described using the supercell cal- 
culation (solid lines) or using the locking technique (dashed 
lines), i.e., using perfect Pt monowires for the leads. 

B. Transmission: scalar-relativistic case 

With the aid of the Pt monowire DFT calculations and 
the construction of WFs and the Hamiltonian from the 
hopping matrix elements it is now possible to calculate 
the conductance based on the Green's function method. 
We start by considering the quality of the locking tech- 
nique. For this purpose, we compare the results for a 
rather small 6 atom supercell calculation for the scat- 
tering region with a single elongated bond of A = 0.72 
bohr and a calculation performed in a 12 atom supercell. 
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The quantum conductance obtained for both cases with- 
out applying the locking-technique, i.e., constructing the 
semi-infinite leads from the supercell calculation, is sim- 
ilar but differs in key details such as a sharp peak just 
below the Fermi energy (compare Fig. ^a.) and (b)). If 
we replace the Hamiltonian for the leads by the one con- 
structed from the MLWFs of a periodic Pt-monowire the 
result changes as follows: While the conductance based 
on the 12 atom supercell calculation is nearly indepen- 
dent on how the lead was constructed, the result for the 
6 atom calculation improves significantly upon using the 
locking technique and is almost indistinguishable from 
the calculation in the larger 12 atom supercell. This 
demonstrates the applicability and quality of the locking 
technique which allows to save a considerable amount of 
computational effort to calculate the ballistic transport 
properties. 

While the previous test has been performed within the 
nearest-neighbor approximation for the tight-binding like 
Hamiltonian we now determine how accurate the calcu- 
lated transmission function is with respect to the num- 
ber of neighbors included. In Fig. [6]Ja) the transmission 
functions, calculated in first, second, and third nearest- 
neighbor approximation and based on the 12 atom su- 
percell for the scattering region with one stretched bond 
of A = 0.72 bohr are presented. The main effect which 
we observe upon including more neighbors is a widening 
of the energy range in which the conductance is non-zero 
as expected from the comparison of the bandstructure 
obtained in the different approximations, cf. Fig. [3j The 
conductance in the vicinity of the Fermi energy which 
is dominated by the localized d-states is well described 
already using second nearest-neighbors. Using only the 
first nearest-neighbor on the other hand results in an off- 
set of the conductance above the Fermi energy which 
originates from a shift of the upper edge of the A3 band 
as seen in the bandstructure. Therefore, we use at least 
the second nearest-neighbor approximation in the follow- 
ing to construct the tight-binding like Hamiltonian. 

In order to understand which states contribute to the 
transmission we can decompose it with respect to the 
orbital symmetry of the Wannier functions. The individ- 
ual transmission channels can be derived from Eqn. (8j), 
by performing the trace operation only over WFs within 
the same symmetry group. In Fig. [6jb) we see that s-d z 2 
states provide an almost perfectly conducting channel in 
a large energy range. Only far below the Fermi energy 
the value drops below 2e 2 /h and in the vicinity of the 
Fermi energy it rises due to the availability of two Ai 
bands, cf. the bandstructure in Fig. [3j The more local- 
ized d xz — d yz states, on the other hand, possess a much 
smaller transmission and their contribution is localized 
in a small energy window. This effect is even more dra- 
matic for the d xy — d x 2_ y 2 orbitals, which show a very 
small overlap and hopping matrix elements leading to a 
sharp peak in the conductance. 

Finally, we turn to the conductance of the Pt monowire 
as a function of the stretched bond length shown in 
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FIG. 6: (color online) (a) Conductance for a Pt monowire 
with a single bond stretched by A = 0.72 bohr based on 
a 12 atom supercell calculation using the first, second, and 
third nearest-neighbor approximation for the construction of 
the transport Hamiltonian from MLWFs and the locking tech- 
nique to attach the leads, (b) Decomposition of the total con- 
ductance (solid line) for the second nearest-neighbor approx- 
imation into the contributions of the s — d z 2 (Ai) (solid red 



line), the d xz —d yz (A3) (dotted blue line) and the d x 
(A4) (dashed-dotted orange line) channels. 
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Fig. [7] For the conductance of a perfect Pt wire, we 
find the expected step-function shape in which each band 
contributes with one conductance quantum Go per spin 
within its bandwidth. Upon increasing the length of a 
single bond in the wire, the overlap between the Wan- 
nier orbitals across the gap decreases, especially for the 
more localized d-orbitals, and as a result the transmis- 
sion drops dramatically. Accordingly, only the contribu- 
tion from the s-d z 2 states survives at large gaps while 
the sharp peak originating from the d xy — d x 2_ y 2 orbitals 
vanishes above A = 1.82 bohr. Another important re- 
sult of this calculation is that the Hamiltonians obtained 
with MLWFs and FSWFs provide nearly the same re- 
sults, i.e., the radial solutions of the FLAPW potential 
are evidently a reasonable choice as FSWFs trial orbitals. 



C. Transmission: spin-orbit coupling 

For heavy transition metals such as Pt spin-orbit cou- 
pling plays an important role and has a significant im- 
pact on the electronic structure. Evidently, the transport 
properties should be equally affected. A suitable method 
to describe the quantum conductance in such systems 
has to be capable of treating SOC. The effect of SOC 
on the electronic structure, namely the coupling of the 
spin quantum number s = h and angular momentum 
quantum number I = 0, 1, 2, ... to the total angular mo- 



mentum quantum number j = h, 



2' 2 



can be seen in 



Fig. [HJ Compared to the scalar-relativistic calculation, in 
which SOC is neglected, Fig.|3j the bandstructure includ- 
ing SOC changes significantly, Fig. [8] In the chain ge- 
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FIG. 7: (color online) Conductance for a nonmagnetic Pt 
monowire in the scalar-relativistic approximation, i.e., ne- 
glecting spin-orbit coupling, with a single bond stretched by 
A. The second nearest-neighbor approximation has been used 
for the transport Hamiltonian. The WFs and hopping matrix 
elements have been constructed from a 12 atom supercell and 
the leads were described by the locking technique. Curves are 
shown for MLWFs (solid lines) and FSWFs (dotted lines) for 
A = 0.0, 0.34, 0.72, 1.22, 1.82 and 2.52 bohr (from left to 
right). 
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FIG. 8: (color online) Bandstructure of an infinite nonmag- 
netic Pt monowire including spin-orbit coupling, (a) Band- 
structure from the FLAPW calculation (big dots) and us- 
ing the Hamiltonian from FSWFs within the first, second, 
and third nearest-neighbor approximation, (b) Conductance 
based on FSWFs for 1 st (dotted line), 2 nd (dashed line) and 
3 rd (solid line) n.n. approximation. 



ometry, the states arc eigenfunctions to the z-component 
(chain axis) of the total angular momentum and we can 
classify the bands by the absolute value of rrij as shown in 
Fig. [S] (a). Thereby, spin-orbit coupling leads to several 
avoided crossings in the bandstructure, e.g. of a s — d z 2 
and d xz jd y z -band around 3 eV below the Fermi level. 
With respect to the scalar-relativistic bandstructure, we 
also observe a significant shift of the d xy and d x 2_ y 2- 
bands towards the Fermi energy. As this band touches 
the Fermi energy at k = ^ the conductance jumps from 
4 Go in the scalar-relativistic case to a value of 5 Go- This 
finding already demonstrates the importance of SOC for 
quantum transport calculations in such systems. 

The general form of the conductance in presence of 
SOC changes significantly, too, due to the lifted degen- 
eracies of bands with different \ttij [-values (see Fig.[8](b)). 
While the conductance at the Fermi level is enhanced 
upon taking SOC into account, the degeneracy of the d xy 
and d x 2_ y 2 -bands in the SR case leads to a higher con- 
ductance of 6 Go below the Fermi energy. Another key 
difference due to SOC is the larger number of steps which 
appear in the conductance as a result of the anti-crossings 
in the bandstructure, in particular, in the energy range of 
3 eV to 1 eV below the Fermi level. In Fig. [9] we also dis- 
play the evolution of the conductance upon stretching a 
single bond in the Pt monowire. Similar to the SR case, 
we observe a rapid decrease of the conductance due to 
more localized d-orbitals. However, due to the spin-orbit 
split bands there is a more pronounced peak structure in 
the conductance. In particular, we find a sharp peak just 
below the Fermi energy which decays more slowly than 
in the SR calculation where it is located slightly lower in 



energy. Our calculations of the conductance are in good 
agreement with those obtained based on fully relativistic 
ultrasoft pseudopotentials and a scattering approach to 
obtain the conductance^. 




E-E F (eV) 

FIG. 9: (color online) Conductance for a nonmagnetic Pt wire 
with a single stretched bond dp t + A including spin-orbit cou- 
pling calculated within a 6 atom supercell and using locking to 
semi-infinite Pt leads based on the Hamiltonian obtained from 
MLWFs in second nearest-neighbor approximation. From left 
to right: one bond stretched by A = 0.0, 0.34, 0.72, 1.22, 1.82 
and 2.52 bohr. 
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FIG. 10: (color online) (a) Majority and (b) minority band- 
structure for a ferromagnetic Co monowire with dco = 4.15 
bohr calculated within FLAPW (big dots) and FSWFs in 1 st 



(dotted lines), 2 n 
approximation . 



(dashed lines) and 3 r n.n. (solid lines) 



IV. Co MONOWIRES 

Another important aspect in transport through 
nanoscale structures is the effect of spin-polarization and 
magnetic order. Due to the reduced coordination num- 
ber in nanostructures the density of states is enhanced 
and according to the Stoner model the tendency towards 
magnetism increases. The reduced symmetry also results 
in a much larger magneto-crystalline anisotropy energy 
(MAE) as the orbital moments become more significant. 
E.g. freestanding and suspended chains of id- and 5d- 
transition-meta ls bec ome magnetic and show giant val- 
ues of the MAE^inl and the effect of colossal magnetic 
anisotropy has been reported^. Here, we demonstrate 
that our method allows spin-polarized transport calcu- 
lations. We consider a simple model system, i.e., a Co 
monowire with a single stretched bond and allow a par- 
allel and antiparallel alignment of the magnetization on 
the two Co electrodes. We calculate the conductance in 
both configurations and determine the magnetoresistance 
as a function of electrode separation. The calculations in 
the antiferromagnetic configuration of the Co monowire 
can also be compared to calculations by Smogunov et al. 
based on the scattering approach and pseudopotentialiP^. 



A. Magnetoresistance 

Compared to the non-magnetic Pt bandstructure, the 
Co chain exhibits a smaller band width due to more local- 
ized 3d-states, and a large exchange splitting (Fig 
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The exchange splitting leads to a net spin moment in 
the unit cell of 2.13/ib- A good overall accuracy in re- 
producing this bandstructure based on FSWFs can be 
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FIG. 11: (color online) Conductance between two ferromag- 
netic Co monowires separated by a gap, A, in (a) parallel and 
in (b) antiparallel alignment of the magnetization. A super- 
cell of 16 atoms has been used for the scattering region and 
the transport Hamiltonian was constructed based on FSWFs 
in the second n.n. approximation. From left to right: gap 
of A = 0.0, 0.45, 1.05, 1.85 and 2.85 bohr. Upper and lower 
part of the plots show the spin-up and spin-down transmission 
channel, respectively. 



achieved if we go up to third nearest-neighbor hoppings. 
For the d-bands and the s-d z i bands around the Fermi 
energy even the second nearest-neighbor approximation 
is sufficient. From the spin-split bandstructure we ex- 
pect a larger conductance in the parallel magnetization 
alignment due to the overlap between minority bands of 
A3 and A4 symmetry. This notion is confirmed by the 
calculated conductance in the two magnetic configura- 
tions as a function of gap size as shown in Fig. [TTJ At 
the Fermi level, we observe majority and minority spin 
conductances of G ma j = e 2 /h and G m i n = 6e 2 /h, re- 
spectively, for a perfect ferromagnetic Co monowire (see 
Fig 11 (a)). As the central bond is stretched the minority 



conductance drops rapidly because it originates from the 
more localized d XZiVZ - and rf X y,x 2 -y 2 -states. The majority 
conductance, on the other hand, is due to s — d^-states 
and decays much more slowly. 

In the antiparallel alignment (Fig [IT] (b)), the conduc- 
tance is the same in both spin channels. There is only 
a small energy window between 1 eV and 1.7 eV below 
the Fermi energy in which the d XZlVZ - and z 2_ v 2-states 
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FIG. 12: (color online) Conductance at Ef for the antipar- 
allel (red filled circles, solid line) and the parallel (green 
solid squares, solid line) alignment of two ferromagnetic Co 
monowires as a function of separation. For the parallel case 
the decomposition into majority (blue open squares, dashed 
line) and minority (orange open squares, dotted line) contri- 
butions is given. The inset shows the BMR ratio as a function 
of separation. 

overlap and at the Fermi energy, the conductance is dom- 
inated by the s — d z 2 -states. The conductance in the an- 
tiparallcl alignment can be interpreted as an envelope of 
spin-up and spin-down transmission functions calculated 
for the parallel case as an electron can only be transmit- 
ted if there are states of the same symmetry in both spin 
channels. The conductance in this configuration is also 
in good agreement with that reported by Smogunov efc 
al. using a scattering approach and pseudopotentials^. 

Based on the obtained quantum conductance at the 
Fermi level we can calculate the ballistic magnetoresis- 
tance (BMR) upon stretching the central bond. The 
BMR is defined as the difference between the conduc- 
tance in the parallel and antiparallel alignment divided 
by the antiparallel conductance: 



BMR: 



Gp{E F ) — Gap(Ef) 
Gap(Ef) 



x 100% 



(23) 



Fig. [12] displays the evolution of the spin-resolved con- 
ductance as a function of gap size for the two magnetic 
configurations. As noted above, the parallel alignment is 
characterized by a rapidly decreasing minority spin con- 
ductance and a nearly constant majority spin contribu- 
tion. However, the minority spin conductance dominates 
until the end of the bond length range which we con- 
sidered. In the antiparallel alignment, the conductance 
of both spin channels is the same and behaves similar 
to the majority spin channel of the parallel alignment 
as it is due to s — d z 2 -states. From this analysis of the 
channel contribution we can understand the fast drop of 
the BMR found upon stretching (inset of Fig. 12 1 of the 
central bond in the monowire. 



V. SPIN-ORBIT SCATTERING AT 
IMPURITIES 

In the previous sections we applied our quantum trans- 
port code to systems with strong spin-orbit coupling (Pt 
monowires) and high spin-polarization (Co monowires). 
In the following we combine the two effects in order to 
study the scattering at impurities in the presence of spin- 
orbit coupling. We consider two types of model systems. 
We begin with non-magnetic Pt monowire with a single 
Co impurity atom and calculate the dependence of the 
conductance on the magnetization direction of the Co 
atom. An analysis of the orbital decomposed transmis- 
sion function allows us to study the influence of SOC on 
the different channels. We find that band mixing due 
to SOC has a pronounced influence, in particular, on the 
contribution from the band with A^symmetry. As a sec- 
ond system, we consider a ferromagnetic Co monowire 
with a single Pt impurity atom and compute the conduc- 
tance for the two magnetization directions of the Co wire, 
either along the direction of the wire or perpendicular to 
it. From our calculations of the conductance including 
spin-orbit coupling we can also determine the ballistic 
anisotropic magnetoresistance (BAMR), i.e., the differ- 
ence of transmission between a magnetization parallel to 
the current and perpendicular to the current. 

While our systems are idealized they can be seen as 
prototypical for experiments which may be performed for 
example by scanning tunnelin g micro scopy in the contact 
regime 7 or in break junction d 1 ! 10 ! 71 !. Scalar-relativistic 
calculations, i.e., neglecting SOC, in a similar geometry 
for a Ni impurity in a Au monowire have been performed 
before^!. 



A. Magnetic impurity in a non-magnetic wire 

We begin our investigation of spin-orbit scattering at 
an impurity by considering a single Co atom in a Pt 
monowire. This is the simpler of the two systems due 
to the nonmagnetic Pt leads. We have already discussed 
the conductance of Pt monowires with and without spin- 
orbit coupling in section |III| Here, we study the con- 
ductance for different magnetization directions of the Co 
impurity atom in order to calculate the so-called ballistic 
anisotropic magnetoresistance (BAMR), which has been 
predicted based on DFT calculations^ and was experi- 
mentally reported for Co break junctions^. 

Before we discuss the calculated conductance, we fo- 
cus on the magnetic properties of our system. From the 
DFT calculations in the scalar-relativistic case we obtain 
spin moments of 2.46/xb for the Co atom which induces 
Pt spin moment of a magnitude of up to 0.27/!^, oscillat- 
ing in sign as a function of separation from the Co atom. 
A similar behavior was found upon including SOC in the 
calculations for both magnetization directions, with a Co 
spin moment of about 2.49^. Including spin-orbit inter- 
action in the calculations gives rise to finite values of the 
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FIG. 13: (color online) Conductance of a Pt monowire with a single Co impurity in (a) the scalar-relativistic (SR) approximation 
and including spin-orbit coupling (SOC) for a magnetization (b) along the chain axis (z) and (c) perpendicular to it (r). In 
addition to the total conductance (black thick line) each panel shows the transmission for a perfect Pt monowire (dashed-dotted 
line) and orbital-decomposed into the Ai- (red dotted line), A3- (blue dashed line), and A4-band (green solid line) contribution. 
The projection onto the spin-up and spin-down states is given for two different directions of the {/-axis, respectively. Below 
each conductance panel the density of states (DOS) is displayed in the corresponding electronic configuration, i.e., SR or SOC, 
for a perfect Pt monowire, the Pt atom adjacent to the Co impurity, the Co impurity, and a perfect Co monowire. The DOS 
is orbital-decomposed similarly to the transmission. 



orbital moments of the atoms, which play an important 
role in determining the energetically favorable direction 
of the magnetization.^ I n our system, the orbital mo- 
ments of the Co atoms are much larger than those of 
the surrounding Pt atoms, and constitute 0.12/is and 
0.19/1 ^ for the magnetization along the chain axis (z) 
and perpendicular to it (r), respectively. Accordingly}-^ 
this results in an energetical preference of the in-chain 
magnetization direction over the out-of-chain direction, 
with a calculated magneto-crystalline anisotropy energy 
(MAE) of 4.3 meV per magnetic atom. 

We now turn to the calculated conductance presented 



in the three top panels of Fig. 13 for the scalar-relativistic 
case and upon including spin-orbit coupling for the two 
different magnetization directions. For reference the or- 
bitally decomposed conductance and the density of states 



(DOS) of a perfect Pt monowire is given in each of the 
three plots and in the panels below, respectively. As 
a general trend, the introduction of a Co scatterer re- 
sults in a non-perfect matching between the spin-split Co 
3d-states and the more delocalized Pt 5d-states (cf. the 
bandstructures in Figs. [3] and 10). In all three 
clear signature of the exchange-split Co A3-band can be 
observed in the overall conductance, most clearly visible 
in the spin- and orbital decomposition. As expected, the 
A4-bands are shifted towards the Fermi energy upon in- 
cluding spin-orbit coupling, but due to the energetical 
mismatch between the Co and Pt A4-bands in SR and 
for both magnetization directions with SOC, this band 
plays only a minor role in the overall conductance. 

Nevertheless, there is a considerable difference be- 
tween the conductance at the Fermi level in the scalar- 
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relativistic case, Gsr = 1.40 Go, and upon including 
SOC either for z-magnetization, Gti — 2.25 Go, or r- 
magnetization, G± = 2.10 Go, as seen in Fig. [Mj The 
main reason for this large difference between SR and 
SOC conductances can be found in the Ai band of SR 
Pt. In this channel the DOS is reduced compared to the 
SOC cases at the Fermi energy at the Pt n.n. atoms and 
there is a corresponding reduction of the conductance, as 
shown in Fig. [Tl The difference of Gn - G± = 0.15 G 



between the two different magnetization directions can 
be found in the larger minority A 3 -state contribution of 
the parallel aligned axis. Here, the SR and the parallel 
SOC case behave similarly. The DOS for A3 majority 
states is small at the Fermi level, the majority state con- 
ductance is reduced in comparison to the minority state 
contribution, as a result of the exchange splitting of the 
Co scatterer. 

Interestingly this is not the case for the r direction of 
the magnetization, for which majority and minority chan- 
nels contribute equally to the total conductance. This 
effect also occurs for the A 3 minority channel between 
—2.8 eV and —3.9 eV as well as for the A4 conductance 
just below the Fermi energy. While the very sharp spin- 
up A4-peak in the SR transmission at —0.7 eV can be 
traced back to a small spin-up A4 peak in the DOS of 
the central Co atom at this energy, this is not the case for 
the mentioned regions in case of the r-magnetization, for 
which no majority A3 and A4 states are present at the 
scatterer. The origin of this effect is the broken cylin- 
drical symmetry when the magnetization points out of 
chain. This broken symmetry allows for a hybridization 
between Ai and A3 bands with j — h, as well as be- 
tween A3 and A4 bands with j = |. As a result an 
incident electron of j = \ (|) can be transmitted into 
a state with j = \ (!) of different orbital character and 
spin. This channel for scattering is less effective than 
the spin-conserving scattering for the in-chain magneti- 
zation, resulting in a larger conductance in this case. 

The changes in the ballistic conductance due to 
ballistic spin-scattering are important for the ballistic 
anisotropic magnetoresistance (BAMR). The BAMR is 
defined analogously to the anisotropic magnetoresistance 
as: 



BAMR = — L — - x 100%, 
G 1 



(24) 



where Gn and Gj_ are the conductances for the magne- 
tization along the wire axis and perpendicular to it, re- 
spectively^. The difference of 0.15 Go at the Fermi level 
in favor of the parallel quantization axis due to ballis- 
tic spin-scattering results in a small BAMR of the order 
of 7%, see inset of Fig. 14 A small shift between the 



A4 contributions due to a small spin-splitting of those 
bands for the Pt atom neighboring the Co scatterer (cf. 
Fig. 13 ) results in an oscillatory behaviour of the BAMR 
when the energy is varied from —0.05 eV to —0.2 eV, 
with BAMR ranging from -20% to 25%. 




E-E (eV) 



FIG. 14: (color online) Conductance around the Fermi level 
for a Pt wire with a single Co impurity atom without spin- 
orbit coupling (dotted blue line) and including SOC for an 
in-chain (dashed red line) and an out-of-chain (solid black 
line) magnetization direction. The inset shows the ballistic 
anisotropic magnetoresistance (BAMR) as defined by Eqn. 24 



B. Non-magnetic impurity in ferromagnetic wire 

In the previous example, we have seen that the trans- 
mission can be affected by ballistic spin-scattering, lead- 
ing to a small BAMR below the Fermi energy and BAMR 
oscillations due to a shift in the A4 orbitals of the Pt atom 
next to the Co scatterer. In this section, we consider a 
non-magnetic scatterer, a Pt atom, in a ferromagnetic 
Co monowire. We find that this situation leads to an en- 
hanced BAMR close to the Fermi level, which is crossed 
by the A4-band. In this case we do not expect strong bal- 
listic spin-scattering because of the magnetic leads, since 
large exchange splitting prohibits scattering between the 
states with opposite spin. 

First we consider the junction in the scalar-relativistic 
approximation in order to understand the main impact 
of the Pt scatterer on the conductance. While Co atoms 
in the leads carry a magnetic moment of 2.13/xb, the 
Co atoms in the vicinity of the Pt atom have moments 
in the range of 2.15— 2.20/is, and the Pt atom itself is 
spin-polarized with a considerable moment of 0.36/is. 
As can be seen in the orbitally decomposed conductance, 
Fig. 15 |a-c), the reduction of the transmission due to 
the Pt impurity atom is relatively small compared to the 
perfect ferromagnetic Co monowire. We can understand 
this general behavior from the fact that the Pt 5d-bands 
possess a broader bandwidth and thereby allow trans- 
mission in the entire regime of the spin-polarized Co 3d 
bands (cf. Fig. [3] and Fig. [To]). 

In the s—d z 2 -channel, the reduction of the transmission 
is similarly small for the majority and minority spin con- 
tributions due to the energetic alignment of the spin-split 
states of the Co wire with the states of the Pt impurity. 
In the majority spin channel, a significant reduction of 
transmission only occurs in a region from Ep — 2.1 eV to 
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FIG. 15: (color online) (a-c) Orbital decomposition of transmission through a ferromagnetic Co wire with a single Pt impurity 
(solid lines) and a perfect Co wire (dashed lines) for the Ai, A3, and A4 channels, respectively, for spin-up (black, upper part) 
and spin-down (red, lower part), (d-f) Density of states of the Co leads orbitally decomposed for spin-up (dashed black lines, 
upper part) and spin-down (red dashed lines, lower part) states, (g-i) Orbitally decomposed DOS of the Pt impurity atom, for 
spin-up (black, upper part) and spin-down (red, lower part). 



Ep — 0.9 eV where the perfect conductance amounts to 
Go . In the spin- and orbital-decomposed density of states 
(DOS), Fig. 15 'd) and (g), we also find two resonances at 
the Pt impurity located at 2.8 eV, and 2.3 eV below the 
Fermi energy in the majority and minority spin channel, 
respectively. In the conductance, we observe a Fano-type 
line shape due to the coupling of the Ai-band to these 
resonances. 

The conductance from the A3-bands displays only a 
reduction at the bottom and top of the band in both 
spin channels as the on-site energies of Co and Pt d xz>yz - 
states are close in energy. The density of states of the Pt 
'h), shows that the d XZjVZ states are spin- 
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atom, Fig. 

split, carry a significant part of the Pt moment, and align 
well with the A3-bands in the Co monowire resulting in 
an efficient transport channel. The most severe change in 
the conductance upon introducing a Pt impurity occurs 
in the A 4 -band. Here, we observe a large decrease due 
to scattering at the Pt impurity. For both, the d xz>yz 
and d xy _ x 2_ y 2 channels, bound states on the Pt atom can 
be found due to the lower on-site potential at the Pt 
site. For the A 3 -symmetry there are such states at —4.1 
eV for the majority band and at —3.3 eV for both spin 



channels, which do not contribute to the conduction as 
they are below the A3-band of the Co leads. For the A4 
symmetry there are majority states around —2.5 eV and 
a paired state at —2 eV with respect to the Fermi level, 
not contributing to the majority channel transmission. 

A4 electrons are only transmitted in the small overlap 
region around —2.1 eV for majority and around Fermi 
level for minority states, where a very narrow band is 
formed in both cases. The shape of the transmission 
function follows the two-peak (majority band) and three- 
peak (minority band) shape of the DOS of the central Pt 
atom. 

Now we turn to the effect of spin-orbit coupling on the 
magnetic and transport properties of the Co-Pt-Co junc- 
tion. For the perfect Co monowire we found a magneto- 
crystalline anisotropy energy (MAE), i.e., the difference 
in energy for the magnetization in the chain axis and 
perpendicular to it, of 0.8 meV per magnetic atom in 
favor of an out-of-chain magnetization and orbital mo- 
ments of 0.17/is for out-of-chain and 0.22/is for in-chain 
direction. Upon introducing the Pt atom, this value is 
reduced to 0.5 meV per magnetic atom, which is con- 
sistent with our observation in the previous section for 
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FIG. 16: (color online) Conductance including SOC for a 
Co monowire with a Pt impurity atom for (a) magnetization 
along the wire axis and (b) perpendicular to the axis. The 
decomposition of transmission into Ai (dotted red lines), A3 
(dashed blue lines) and A4 (solid green lines) channels for 
majority spin (positive y-axis) and minority spins (negative 
y-axis) shows the presence of a A4 minority band channel at 
the Fermi level. Black dashed-dotted lines display the trans- 
mission of the perfect infinite Co leads. 



a Pt-Co-Pt junction favoring the in-chain direction. The 
magnetic moment of the Pt atom is 0.36/zs for both mag- 
netization directions and we find similar orbital moments 
of 0.09 fiB (out-of-chain) and 0.10/zb (in-chain). Charac- 
teristically, as in the case of the Co leads, the orbital 
moments of the Co atoms adjacent to the Pt impurity 
are significantly larger for the in-chain direction (reach- 
ing as much as 0.53/xs for the nearest Co atom), than for 
the out-of-chain direction (at most 0.2/is). This means 
that the out-of-chain easy magnetization axis in our scat- 
tering region is mainly due to the Co atoms. 

For the transport properties including SOC, the bands 
with A 3 and A4 symmetry are essential. Depending on 
the quantization axis defined by the magnetization direc- 
tion, the degeneracy of these bands is lifted. In contrast 
to the Pt-Co-Pt system, the Co electrodes are ferromag- 
netic and therefore the splitting in the step-like conduc- 
tance in the perfect Co wires changes upon switching the 
quantization axis from along the chain axis to perpendic- 
ular to it. 

As can be seen from Fig. |16[ changing the magnetiza- 
tion direction in a perfect infinite Co chain leads to a re- 
duction of the transmission from 3.5 Go (along the chain) 
to 1.5 Go (perpendicular to the chain) in a very small 
energy window around the Fermi energy, which results 
in a huge value of the ballistic anisotropic magnetoresis- 
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FIG. 17: (color online) Conductance around the Fermi level 
for a Co monowire with a Pt impurity atom in the scalar- 
relativistic approximation (blue dotted line) and including 
SOC for a magnetization along the chain axis (dashed red 
line) and perpendicular to the axis of the wire (solid black 
line). The inset shows the BAMR as defined by Eqn. ( 17 1 for 



the Co monowire with a Pt impurity (solid black line) and for 
a perfect infinite Co monowire (dashed-dotted green line). 

tance of 133%P^ In a realistic situation, however, such 
values of the anisotropic magnetoresistance can be hardly 
achieved, owing to the destruction of perfect conducting 
channels by imperfections, impurities and disorder. 

In the case of a Co chain with a Pt impurity, similarly 
to the scalar-relativistic case, we observe a reduction by 
roughly a factor of two in the overall conductance over 
the entire energy range, due to the less efficient coupling 
between the Co wire and the Pt impurity compared to an 
perfect Co wire, especially for the A 3 - and A4-orbitals. 
At the Fermi energy, we find majority and minority spin 
contributions from the A 1 -band of about 0.5 and 1.0 Go 
for both magnetic directions. Only the minority states of 
the other two orbital symmetries are present due to the 
exchange splitting. The minority A3 band contributes 
almost 1.0 Go for the in-chain magnetization, while it 
reveals a large dip at Ep for the out-of-chain magnetiza- 
tion. Accordingly, the A4-band conductance also changes 
significantly upon switching the magnetization direction, 
owing to the changes in the details of hybridization be- 
tween A3 and A4 states when the direction of the mag- 
netization is changed, see Fig. 16 'b) and Fig. 13 (cf. DOS 



of the Co monowires for the two different magnetization 
directions). These changes in the energetic structure of 
A 3 and A4 states lead to a large difference between the 
in-chain and out-of-chain conductances, also visible for 



the pure Co chain in the Fig. 16 



In Fig. [TTJ the conductance is displayed in a small en- 
ergy window around the Fermi energy for the two dif- 
ferent magnetization directions. It is apparent that the 
changes arise due to the modifications of the A4-band 
conductance between Fermi level and —0.05 eV, and A 3 - 
band conductance around Ep and —0.15 eV which are 
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subject to different band mixing from spin-orbit coupling. 
As a result of the fine structure of the A4 and A3 con- 
ductances (sec Fig. [l6^b)), the BAMR which we obtain, 
shown in the inset of Fig. |17[ displays a strong varia- 
tion with energy. Compared to the BAMR of a perfect 
Co MW of 133% a Pt scatterer reduces this effect to 
80 — 100%, which is still considerably high. An enhanced 
BAMR can be found for the second peak below the Fermi 
Energy, where a A4 conduction peak for the in-chain di- 
rection in coincidence with a A3 conduction depletion 
result in a BAMR increase from 40% for the perfect Co 
MW to 60 — 100% when a Pt scatterer is introduced. 



VI. SUMMARY 

We have implemented the Landauer-Biittiker method 
to calculate the ballistic electron transport through one- 
dimensional nanoscale junctions based on density func- 
tional theory calculations within the full-potential lin- 
earized augmented (FLAPW) method. In order to ap- 
ply the efficient Green's function method to calculate the 
conductance we have mapped the extended Bloch states 
obtained from the FLAPW method to the minimal basis 
set of localized Wannier functions and constructed the 
Hamiltonian for the open system. With our approach 
it is feasible to calculate ballistic transport through one- 
dimensional nanoscale systems including magnetism and 
spin-orbit coupling with the accuracy and flexibility of 
the FLAPW method. 

We apply our method to calculate the conductance of 
non-magnetic Pt monowires with a single stretched bond, 
including spin-orbit coupling. Already this simple exam- 
ple shows the key impact of SOC for systems contain- 
ing heavy transition metals. As a second example, we 
considered a Co monowire and studied the magnetore- 
sistance upon stretching the wire at a single bond. The 
decomposition of the transmission into the channels of 
different orbital symmetry shows the dominant contribu- 
tion of s- and d^-states as one moves from the contact 
to the tunnel regime. Finally, we studied the effect of 
spin-orbit scattering at an impurity atom in a monowire. 
We considered two model cases: (i) a magnetic atom in 
a non-magnetic wire, Co in a Pt monowire, and (ii) a 
non-magnetic heavy element in a ferromagnetic wire, Pt 
in a Co monowire. We observed for both cases a distinct 
dependence of the conductance on the magnetization di- 
rection with respect to the wire axis. 

We found for a Co impurity in a Pt chain, that due 
to the broken cylindrical symmetry for an out-of-chain 
magnetization direction the hybridization between states 
of different angular character and spin but with identi- 
cal quantum number j leads to scattering processes that 
do not conserve spin. Those ballistic spin-scattering pro- 
cesses are resulting into a BAMR of 7%. The relatively 
moderate values are caused by the large background con- 
ductance from bands originating from s—d z 2- and d xz ^ yz - 
states which are not modified much upon switching the 



magnetization. On the other hand, for a Pt impurity 
in a Co chain we find that the presence of an impu- 
rity, although reducing somewhat the BAMR of the pure 
Co chain, still leads to values of BAMR of about 100%, 
which originates from hybridization between the A3 and 
A4 states moderated via SOC by the direction of the 
magnetization. 
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Appendix A: Computational details 
1. Pt monowires 

Non-magnetic (NM) 6 and 12 atom supercell calcula- 
tions with an interatomic distance of dp t = 4.48 bohr 
and the central bond stretched by A = 0.0, 0.34, 0.72, 
1.22, 1.82 and 2.52 bohr. We applied the generalized gra- 
dient approximation (GGA) to the exchange-correlation 
potential^. For calculations in the scalar-relativistic 
(SR) approximation, the irreducible part of the ID Bril- 
louin zone (BZ) was sampled by 6—10 fc-points depending 
on the size of the supercell. For the 6 atom supercell, we 
also performed calculations including spin-orbit coupling 
in second variation. For calculations with SOC the whole 
ID BZ was sampled by 24 fc-points. In all calculations, 
G max was chosen to be 3.7 bohr -1 , which corresponds to 
approximately 200 basis functions per atom. The diame- 
ter of the cylindrical vacuum, -D vac , and the value of the 
in-plane auxiliary lattice constant, were set to 5.0 

and 7.3 bohr, respectively. 

For the conductance calculations we applied the lock- 
ing technique to a perfect monowire to describe the semi- 
infinite leads (see Sec. II D I. In the SR approximation 
FSWFs and MLWFs were generated on a mesh of 16 
fc-points in the whole BZ starting from one 4s- and 5 
3c?-orbitals per atom in the supercell, based on solutions 
of the radial equation of the first-principles potential as 
trial functions. In the calculations including SOC, ML- 
WFs were generated on a 24 fc-point mesh in the whole 
BZ based on 2 radial 4s- and 10 radial 3d-orbitals per 
atom, based on solutions of the radial equation of the 
first-principles potential as trial functions, due to the cou- 
pled spin channels. The energy bands were disentangled 
using the procedure described in Ref. [55] • For the SR 
calculations, the lowest 80 eigenstates are needed for 72 
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WFs for the 12 atom supercell and the lowest 44 eigen- 
values per fc-point for 36 WFs for the 6 atom supercell 
calculations. With SOC the lowest 80 eigenstates per 
fc-point for 72 WFs were used. 



2. Co monowires 

Calculations with a lattice constant of dc = 4.15 bohr 
and a central stretched bond with stretching A = 0.0, 
0.45, 1.05, 1.85 and 2.85 bohr. Two collinear magnetic 
configurations of the Co monowire are considered, paral- 
lel or antiparallel alignment of the Co spins on the left 
and on the right sides of the gap, described by performing 
two calculations: A 8 atom supercell constructed from 
two 4-spin blocks separated by a gap and aligned in par- 
allel (up), while in order to mimic the antiparallel align- 
ment, we considered 16 atoms in the supercell with 4-spin 
(up), 8-spin (down) and 4-spin (up) blocks, separated by 
two gaps with the spins antiparallel to each other at each 
side of the gap. 

The perfect lead ferromagnetic Co monowire was cal- 
culated with 24 fc-points in the whole BZ, using the G max 
of 4.1 bohr -1 (sa 220 basis functions per atom). For both 
8 and 16 atom supercell calculations the irreducible part 
of the ID Brillouin zone (BZ) was sampled by 8 fc-points 
and G max was chosen to be 3.7 bohr -1 resulting in ap- 
proximately 210 basis functions per atom. The vacuum 
parameters Z? V ac and D constituted 4.3 and 6.6 bohr, 
respectively, in all cases. The exchange-correlation po- 
tential was treated within the GGA^. For all quantum 
conductance calculations the locking technique (see Sec. 
II D I to a perfect FM Co monowire was used. As trial 



orbitals for the FSWFs 6 s- and d-orbitals per atom and 
spin in the supercell were used, based on solutions of the 
radial equation of the first-principles potential. For the 
disentanglement procedure^ the lowest 58 (110) eigen- 
states per fc-point were used to obtain the 48 (96) WFs 
in the 8 (16) atom supercell calculation. 



3. Scattering on impurities 

A 9 atom supercell were used for the scattering re- 
gion consisting of one impurity atom (Pt or Co) and 4 
monowire atoms (Co or Pt) on both sides. The inter- 



atomic distance was chosen as dc — 4.15 bohr for the 
Co monowire with a Pt impurity and as c?pt = 4.48 bohr 
for the Pt monowire with a Co impurity. The exchange- 
correlation potential was treated within the GGA^ and 
SOC was included in second variation. All calculations 
were performed in the scalar-relativistic (SR) approxima- 
tion and for two different directions of the magnetization 
with SOC, along the chain axis and perpendicular to it. 
The ID Brillouin zone (BZ) was sampled by 16 fc-points 
and G max was set to 3.9 bohr -1 resulting in approxi- 
mately 175 (190) basis functions per atom for the Co 
(Pt) monowire with a Pt (Co) impurity. For the case of 
an isolated Pt impurity, the leads were described by a 
Co monowire in a 3 atom unit cell in either the scalar- 
relativistic (SR) approximation or including SOC for the 
magnetization direction along the wire axis or perpen- 
dicular to it. The BZ was sampled by 24 fc-points and 
G max was set to 4.1 bohr -1 , resulting in approximately 
210 basis functions per atom. For Pt monowire with a 
Co impurity, the lead's electronic structure was obtained 
from calculations of perfect Pt monowires. The vacuum 
parameters for all cases constituted 4.3 and 6.6 bohr for 
D vac and D, respectively. 

For all quantum conductance calculations the locking 
technique wa used and the third nearest-neighbor ap- 
proximation were employed. In the SR case FSWFs were 
generated on a 16 fc-point mesh in the whole 1D-BZ with 
1 s- and 5 d-orbitals per atom and spin, based on solu- 
tions of the radial equation of the first-principles poten- 
tial. For the disentanglement procedure^ the lowest 64 
(62) eigenvalues per fc-point for 54 (54) WFs for Pt (Co) 
impurities in Co (Pt) monowires were considered. The 
Pt and Co lead WFs were constructed as described in 
Sec.EDand Sec. Klin this case. With SOC the FSWFs 
were generated on a 16 fc-point mesh in the whole 1D-BZ 
with 2 s- and 10 d-orbitals per atom, based on solutions 
of the radial equation of the first-principles potential. For 
disentanglement!^ the lowest 116 eigenstates per fc-point 
for 108 WFs were used. The WFs for the semi-infinite 
Co leads were generated on a 24 fc-point mesh with the 
same trial functions as those used for the atoms inside 
the scattering region, while for disentanglement the low- 
est 26 eigenvalues per fc-point for 18 WFs per spin (SR) 
and the lowest 44 eigenvalues per fc-point for 36 WFs 
(SOC) were used. The Pt lead WFs were constructed as 
described in Sec. IA 11 
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